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Abstract 



It has recently been shown that there are efficient algorithms for quantum 
computers to solve certain problems, such as prime factorization, which are in- 
tractable to date on classical computers. The chances for practical implemen- 
tation, however, are limited by decoherence, in which the effect of an external 
environment causes random errors in the quantum calculation. To combat this 
problem, quantum error correction schemes have been proposed, in which a sin- 
gle quantum bit (qubit) is "encoded" as a state of some larger number of qubits, 
chosen to resist particular types of errors. Most such schemes are vulnerable, 
however, to errors in the encoding and decoding itself. We examine two such 
schemes, in which a single qubit is encoded in a state of n qubits while subject 
to dephasing or to arbitrary isotropic noise. Using both analytical and numerical 
calculations, we argue that error correction remains beneficial in the presence of 
weak noise, and that there is an optimal time between error correction steps, 
determined by the strength of the interaction with the environment and the 
parameters set by the encoding. 



1 Introduction 



Soon after the discovery of fast quantum algorithms for factorization jT|, it was realized 
that the efficiency of quantum computers depends crucially on the control of errors 
during a computation. This is not surprising in itself, since classical computers also 
require an active monitoring of errors to operate properly. However, the dissipative 
techniques used in classical error correction destroy the superpositions necessary for 
quantum computation. This problem stimulated an important effort in the direction of 
quantum error correction. In the last year or so, after the initial discovery of quantum 
error correction codes by Shor and Steane || |J, significant progress has been 
made in the development and understanding of these codes. Much attention has been 
devoted to the construction of codes using a variety of different techniques to convert 
classical codes into quantum codes || |6| and providing a mathematical description of 
large families of these codes || . Minimal codes, that correct only one or a few errors, 
were also derived [|, [UJ . Most of this work addresses the issue of how most efficiently to 



preserve a quantum state in a noisy environment given that encoding and decoding can 
be done in an error-free way. Only a few recent papers have addressed the possibility 
of encoding and decoding in the presence of noise |ll|, |12| ; these fault-tolerant schemes 
are relatively complicated and involve many more qubits than the earlier simple codes. 
It is therefore unlikely that we will see an experimental implementation of these more 
elaborate proposals in the near term. On the other hand, the issue of errors arising 
during encoding and decoding has only been partially investigated in the simplest error 
correcting codes proposed so far jl3||. These are the codes that could be implemented 



in a near-future quantum computer. The aim of this work is to bridge this gap, and 
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provide numerical as well as algebraic evidence that for certain regimes of noise, error 
correction is worthwhile even when noise is present during the encoding and decoding 
steps, as will be the case in any real experiment. 

This paper is organized as follows. Section |2] reviews the effect of unwanted 
environmental coupling on a quantum computer, its description in terms of master 
equations, and the fundamental operation of error correcting codes. Section ||| presents 
an analytical discussion of the effect of errors on an evolution consisting of encoding, 
free evolution of an encoded qubit, and decoding, with noise present at all stages. The 
concept of a quantum trajectory is described in section f|, where numerical simulation 
algorithms are presented using two unravelings of the master equation (quantum jumps 
and quantum state diffusion) as well as direct numerical solution of the master equation. 
The numerical methods are found to be in reasonable agreement with the analytical 
model. 



2 Noise and error correction 
2.1 Decoherence 

A quantum system in complete isolation evolves according to the Schrodinger equation 



where \ip) is the state of the system (in this case a quantum computer) and H is the 
Hamiltonian (in this case representing the action of the quantum "gates;" in general, 
H will be time-dependent). This evolution is unitary. 

Unfortunately, the approximation of a system being isolated is only good for 
microscopic noninteracting systems. As a system becomes larger and more complicated, 
the effects of the environment become more important. 

Consider the example of a single qubit interacting with an external environment. 
The state of the qubit is described by a vector in a two-dimensional Hilbert space. A 
convenient basis is the canonical basis B = {|0), |1)}. Suppose that the qubit is initially 
in a superposition state a|0) + /3|1) and the environment in some unknown state \A). 
As the system and environment interact, the initial product state (a|0) + /5|1)) <8> \A) 
can evolve into an entangled state a|0) <S> \Bq) + /3\1) <S> \B\), where the environment has 
become correlated with the state of the system (more realistic models of coupling with 
the environment can be found in Ref 0]). The system can no longer be described by 
a state on its own. Normally, an environment is very complicated, containing many 
degrees of freedom, so it is likely that \B ) and \Bi) will be orthogonal (or very nearly 
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so). Thus, if we trace out the environment degrees of freedom, an ensemble of our 
systems of interest is left in a mixture, described by a reduced density matrix 



p = i«no><oi + mi><i|. (2) 

In effect, the environment has measured the value of the qubit, and the superposition 
has been destroyed. In this case, the evolution of the reduced system is no longer 
unitary; and algorithms which depend on the unitarity of the evolution, such as the 
Shor algorithm, will no longer function. 

The general effects of the environment on a quantum system can be very compli- 
cated and difficult to describe. However, a useful approximation is to assume that the 
effects of the environment are Markovian, or local in time. In this case, it is possible to 
describe the evolution of the reduced density matrix by a master equation of Lindblad 
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J t P = ~[H,p}+T, (LjpL) - \L]L 3 p - \pL]L^ , (3) 

where H is the system Hamiltonian and the Lj are the Lindblad operators representing 
the interaction with the environment. This Markovian approximation is generally very 
good when the environment is large compared with the system and the interaction 
between them is fairly weak. It might fail, however, for some realizations of quantum 
computers. 

What kinds of Lindblad operators typically occur in (0)? This depends on the 
physics of the system, but certain operators are common in quantum optical and atomic 
physics models. One normal effect of environmental interaction is dissipation, as in 
spontaneous emission. If the qubit state |1) represents an excited state, there will be 
a Lindblad operator proportional to the lowering operator, of the form L = ^/K^^-, 
(<r_|l) = |0) and <x_|0) = 0). The qubit will tend to the ground state in the long term, 
regardless of its initial state. If the environment has a non-zero temperature, there is 
the possibility of thermal excitations as well, represented by another Lindblad operator 
proportional to the raising operator <t + = a_ . 

Even if the rate of dissipation is small enough to be negligible, the environment 
can still act to destroy superpositions. As we saw in (TJ), correlations which develop 
with the environment can randomly dephase the basis states |0) and |1). This process 
is represented by a Lindblad operator proportional to the z Pauli matrix, L = ^/~ko- z , 
(£*|0) = -|0)and£,|l) = |l)). 

The most general interaction will reduce the qubit ensemble density operator to 
the one at the center of the Bloch sphere. The individual (pure) states of the members 
of the ensemble can be viewed as moving randomly on the surface of the sphere. This 
effect is represented by isotropic noise, with three Lindblad operators proportional to 
a x ,a y and a z , used in studying the depolarizing channel [0, jlQfl . The exact choice of 
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model is determined by the physics of the quantum computer and its environmental 
interactions. 



2.2 Error correcting codes 

For a qubit the most general form of single-qubit error induced by the environment 
can be written as 

Wl - (4) 
where E is an arbitrary operator which can be decomposed as 

E = eil + e x a x + e y a y + e z a Z} (5) 

where 1 is the identity and the <7j are the Pauli matrices. In the simplest case (when 
one wants to protect a single qubit), error correcting codes consist in encoding the 
basis states |0) and |1) of a qubit in well chosen states of several qubits: 

where |C ) and \C\) belong to the extended Hilbert space of several qubits. Numerous 



techniques for constructing these codes have appeared in the recent literature [|T6 
These error correcting techniques commonly assume that the encoding step (i.e. the 
operation by which a single qubit in state a|0) + (3\1) is entangled with additional 
qubits to form the state a\C ) + (3\Ci)) as well as the decoding and correcting steps are 
done in a noiseless environment. The issue of noisy encoding and decoding has been 
addressed little outside the context of fault-tolerant techniques, which require many 
more qubits ]TI, I2| . In this work, we will focus on earlier and more compact codes, 



and analyze the issue of noisy encoding and decoding. In the next two sections we 
review the two codes that we have analyzed. 



2.3 Dephasing noise 

If one seeks to protect a qubit against a dephasing noise (in Eq. |5] this is equivalent 
to setting e y = e z = 0), it can be shown that the smallest possible code requires 
three qubits to encode the states |0) and |1) of the initial qubit. This carefully chosen 
superposition was first proposed by Shor p. We use here an equivalent version found 
in Ref . |5|, [T^J , in which 



|0> 
|1> 



T7fl , in which 

|C ) = |000) + |001) + |010) + |011) + |100) + 1 101> + |110) + |111) 
\d) = jooo) - |001) - |010) + |011) - 1 100) + |101) + |110) - 1 111> 



(7) 



(normalization factors have been omitted). We use the networks for encoding and for 
decoding/correcting shown in Fig. [I]. Initially, the first qubit is in the state a|0) + 
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This is the state we wish to protect. The second and third qubits are in the state |0). 
The result of the encoding network is the three-qubit state a\Co) + (3\C\). 

In these figures we choose the various quantum gates in such way that they can 
be described by simple Hamiltonians. The gates A correspond to the unitary operation 




effected by the Hamiltonian Ha = f {~^m(px — &z) + 1) acting for one unit of time. 
Similarly the gate U y corresponds to the unitary operation 

(9) 




generated by the Hamiltonian H y = —\o y acting for one unit of time. 

Please note that these matrices are represented in the basis {|0), |1)}, as is the 
convention in the quantum computation literature. Unfortunately, this is precisely the 
opposite of the usual convention for the Pauli matrices. For this paper we have retained 
the quantum computation basis, as we do not present the Pauli matrices explicitly, but 
this notational conflict should be resolved. 

The two-bit gates of the network correspond to controlled phase shifts. These are 
represented in the canonical basis B = {|00), 1 01) , |10), |11)} by the unitary operator 
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(10) 



generated by the Hamiltonian H = 7rP|iw (g> P\i)j, where i and j designate the two 
qubits on which the gate acts and P\i) = |(1 + a z ) is the projector on state |1). In 
this case a state picks up a phase tc iff both qubits are in state |1). Variants on 
this gate can be obtained by replacing either or both of the projection operators with 
Pio) in the definition of the Hamiltonian. 

The decoding network is just the reverse of the encoding network. After complet- 
ing the sequence of gates, qubits 2 and 3 can be measured to identify the error, which 
is followed by an adequate correction of the first qubit [|, O . 



2.4 Arbitrary noise 

In the previous section we have shown how to protect a single qubit against dephasing 
noise (i.e., noise generated by a single Lindblad operator proportional to a z ). If one 
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Figure 1 : Encoding and decoding networks from [[TtJ . The graphical conventions are similar 
to those in Ref. ||i8fl . The gate Ua represents the operation |0) — ► |0) + |1) and |1) — > |0) — |1) 
(cf. Eq. ^ in the text). The gate U y is a it rotation of the qubit along the y axis (Eq. ||). The 
two-bit gates denoted by a line and two black dots are "control-phase" gates. They change 
by 7r the phase of a quantum state only when both qubits are in state |1) (cf. Eq. [lO| ) . 

seeks to protect a single qubit against an arbitrary error of the form fl5|), such as 
isotropic noise, then five qubits are necessary to encode a state. Different versions of 
these codes (all equivalent) have been proposed ||, We choose to implement an 
equivalent version of the code given by ||: 

\C )= |6 1 >|00)-|6 3 )|H) + |&5)|01> + |67)|10) (u) 
|Ci) = - |6 2 )|11)-|6 4 )|00)-|6 6 )|10) + |6 8 )|01), [ii> 



where |6i) = (|000) ± |111)), |6s) = (|010) ± |101>), \bs) = (|001) ± 1 110>), \br) = 

2 4 6 8 

( 1 1 1 ) ± 1 100) ) . The implementation of this code is done in a way similar to the three 
bit case. In the first stage, a qubit in state |^ n ) = a|0) + /3|1) is entangled with four 
additional qubits (initially in state |0)) to produce the state a\Co) + (3\C\). This is 
done through the network of Fig. 0. Similar gates as in the previous section are used. 
Note that some of the control-phase gates change the phase when a qubit is in state 
|0) rather than in state |1), unlike the three qubit code. 

The decoding network is also given by Fig. [| with the gate operations performed 
in the reverse order (this is possible because each gate appearing in the network is self- 
adjoint). After decoding, qubits 2 to 5 are measured and, depending on the outcome, 
an appropriate correction is applied to return the first qubit to the correct state. (For 
a complete description of this code, see Ref. ||.) 



3 Analytic considerations 

Given that error correction can be implemented through encoding a one-qubit state into 
an n-qubit state, it is instructive to consider some simple analytic conditions for the case 
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Figure 2: Encoding network for the five bit error correcting code from Ref. []|]. The network 
has been adapted in order to provide a more natural implementation of the gates in term 
of Hamiltonians. The decoding network is identical to the encoding one read backwards. 
Dashed line indicate units of time, i.e., all the gates between two dashed line can be effected 
in one unit of time. The total encoding (or decoding) time is therefore 10 units. The gate 
that appear at time 2 is a generalization of the control-phase gate; it imparts a phase factor 
of 7r to the quantum state only if all qubits denoted by a black dot are in state |1). It will 
effect the operation |1, ig, 1, 1, 15) — ► — 11, £2, 1, l,is) and leave all others states untouched. 
The gate at time 3 is a variation: the qubits indicated by a white dot act as controls when 
they are in state |0) rather than |1). 
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of imperfect encoding, decoding and correction. These conditions will complement our 
numerical results and indicate the parameter ranges for which correction is potentially 
useful. 

Subsections pTT] - pT3| discuss a simple approach. To model the imperfect operation 
of the code, we assume that the decoherence which acts during the encoding and decod- 



ing does not get corrected at all. Subsection examines the validity of modeling the 
influence of the environment as instantaneous errors of the type (|J), and demonstrates 
that this type of treatment can be made consistent with the master equation (|]) to 
low order. 



3.1 Perfect single error correction 

First we need a benchmark. For simple decoherence, the probability that a single qubit 
remains error free for a time T is defined to be 

Psnc = e" K " T • (12) 

The "s" stands for success — it could be that the aim is to successfully store the qubit 
for time T, or to transmit it down a channel where the time taken for this is T — 
and "nc" indicates that no correction procedure is applied to the qubit. If £ is the 
probability of an error and ^<1, then £ = 1 — p snc pa K n T . 

The subscript n on the decoherence rate n n denotes the fact that the number of 
qubits needed for encoding is determined by the type of noise. The relevant examples 
presented in the previous section are n = 3 for phase noise (modeled by L x = \J~kct z ) 
and n = 5 for isotropic noise (modeled by L\ = \fk~a x , L 2 = \l~Kb y and L 3 = \J~k~a z ). 
We define the mismatch between the ensemble at time T and that T = to be 

m nec (T) = 1 - (lpini\p\lpini) , (13) 

where \ipmi) is the initial pure state at time T = 0. In the phase noise example, 
an ensemble of single qubits given by \ipini} = 2 -1 / 2 (|0) + |1)) decoheres so that the 
mismatch m nec is easily shown to be 

m nec (T) = i (l - e- 2KT ) . (14) 

We therefore obtain K3 = 2k, since the exponentially decaying term in the mismatch can 
be identified with the probability that the system remains error free. In the isotropic 
noise case, the same initial ensemble decoheres and exhibits a mismatch of 

m nec (T) = i(l-e- 4ttT ) (15) 
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and so we identify k 5 = 4k. The particular choice of initial state is made so that 
it exhibits sensitivity to phase or to isotropic noise. It is also used in the numerical 
simulations, so we compare like with like. 

Consider now the n-qubit encoding and decoding procedure which is able to cor- 
rect perfectly for a single error in one of the n qubits, but fails if there are two or 
more errors. Using this procedure, the probability of the successful survival of a single 
encoded qubit state for time T is the sum of the zero error and one error probabilities; 

p sc (n) = e~ nKnT + ne - (n - 1)K " T (l - e~ KnT ) . (16) 

(Each qubit is assumed to suffer the same decoherence rate n n .) Clearly, if K n T <C 1, 
perfect error correction is worthwhile because p sc (n) is closer to unity than is p snc . 
(Psnc ~ 1 — K n T, whereas p sc ~ 1 — "^ n ~ 1 - > k^T 2 .) For the case n = 3 the crossover point 
arises when p snc = p sc (3), which yields n 3 T = In 2. For the case n = 5, p snc = p sc (5) 
gives k 5 T = 0.14. Any realistic systems are likely to be well down in n n T from these 
values and so would certainly benefit from perfect error correction. 

3.2 Imperfect error correction 

Consider now the case where the encoding (E) and decoding (D) for the error correction 
procedure take a finite amount of time. Let this be TS, so S is the dimensionless fraction 
of time taken by one full E+D stage. D may be essentially the inverse of E and so each 
may take 5/2, but this is not crucial. The decoherence rate could well differ during 
E and D; we denote it by K,' n . The environment seen by the qubits may be different 
when the encoding and decoding interactions are occurring. The point to note is that 
errors which occur during the E+D stage are unwelcome. We assume that they don't 
get corrected and so contribute directly to the error rate for qubit system. 

If the problem at hand is the storage of a given qubit state for time T, it seems 
reasonable to allow E+D to be part of T, so the encoded ra-qubit state is then kept for 
(1 — 5)T. Alternatively, if the goal is to propagate a qubit state down a channel where 
the time for transmission is T, it would seem to be more realistic to add on T5, so the 
whole process takes (1 + 5)T. This distinction does not appear to be crucial, but we 
examine both cases. 

3.2.1 Storage with imperfect correction 

The probability s that there is no error in this system is the product of the probability 
of all n qubits surviving T5 at a decoherence rate of n' n with the probability of zero or 
one error (which can be corrected) during the time (1 — 5)T at a rate n n . Using (|16|), 
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this gives 



This reduces to p sc (n) given in (|16|) as 5 — > 0. 

To make a simple comparison, let «4 = k„. Equating s sc (n) to p snc {n), the aim is 
to find the crossover value for S. Clearly 5 is about 1/n ; the next order correction in 
n n T gives 

S = - - ^f^nT + ... (18) 

Provided that 6 stays below this value, there should be benefit from error correction 
even though errors may occur during E+D. However, if 5 exceeds this value, the per- 
formance of the procedure is actually worse than doing no correction to a single qubit. 
For the cases n = 3 and n = 5 and provided that K n T <C 1, the bounds on 5 are not 
very constraining. Practical systems would probably have 5 <C 1 and so would operate 
in the regime where imperfect correction is beneficial. 



3.2.2 Transmission with imperfect correction 

The probability t that there is no error in this system is the product of the probability 
of all n qubits surviving T5 at a decoherence rate of n' n with the probability of zero or 
one error (which can be corrected) during the transmission time T at a rate K n . Using 
(|16|), this gives 

t sc ( n ) = e~ n < TS (ne- (n - 1)K " T - (n - l)e~ nKnT ) . (19) 



This also reduces to p sc (n) given in (|i6|) as 5 — >• 



To again make a simple comparison, let n' n = K n . Equating t sc (n) to p S nc{n), the 
aim is to find the crossover value for 5. Once again S is about 1/n ; the next order 
correction in K n T gives 

6= i - ^±K n T + ... (20) 
n 2 

The conclusion for transmission is the same as that for storage. Practical systems with 
n n T <C 1 and 5 <ti 1 will be in the regime where correction is beneficial. 



3.2.3 Single correction optimization 

In the numerical simulations presented in Sect. U the encoding and decoding take a 
set amount of time, rather than a set fraction of T. We therefore define an alternative 
parameterization of the time taken for E+D, setting T5 = A. Our analytic expressions 
can be viewed either in terms of 5 or of A, whichever is most appropriate. 
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Here we just give a simple analytic result, optimizing the time T to achieve the 
most benefit from correction. Assuming that the time taken for E+D is fixed at A, 
what is the optimum T? We find this by maximizing the ratio R of the mismatch 
without correction to the mismatch with correction 

^i^jwM (21) 

1 - s sc {n) 

and similarly for the transmission case. The maximum of R is an indicator of where 
(imperfect) error correction is giving the maximum benefit in comparison to performing 
no correction at all; it is one of the measures we use in our numerical work. To the 
lowest approximation (assuming that K n T <C 1 always), the optimum T is the same 
for both storage and transmission, and is given by 

/ 2A y , , 

Obviously, if the error correction is perfect (effectively taking zero time), we arrive at 
the conclusion that it should be performed as often as possible; T opt = 0. However, for 
cases of practical interest (finite A) this is not so as T opt is then finite. The dependence 
of T opt on k will be compared to our numerical simulations in Fig. |6|. 



3.3 iV-correction procedure 

The basic aim of error correction (within the context of this paper) is to maximize 
the probability of success, storing or transmitting the state as well as possible. The 
parameters T, and 5 are therefore set by the problem at hand. T is set by the 

total length of the transmission channel or the total required storage time. (The latter 
might be the time for which the state "idles" between interactions in a larger quantum 
computation; for such a case the simple error correction procedures discussed in this 
paper would be sufficient to keep it coherent while it idles.) The decoherence is set by 
the environment, n and 5 will be set by the chosen correction scheme and its physical 
realization. However, there is still some freedom. Given all the parameters above, the 
number of correction procedures applied during T can be varied. 



3.3.1 Perfect error correction 



Consider then the problem of optimizing error correction to achieve the greatest prob- 
ability of successful storage or transmission of a qubit state, given the freedom to apply 
an arbitrary number N of E+D procedures during the time T. Assume that these are 
spaced out equally. For the case of perfect error correction, where 5 = so there is 
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no possibility of an error occurring during E+D, it is obviously beneficial to apply as 
many corrections as possible. The probability of success for N applications is 

PNsc(n) = (ne-W*nT/N _ 1)e -™„ t/n^ (23) 

This maximizes for N — > oo, tending to unity independent of the value of n, and is 
consistent with our observation that T opt = when A = 0. Such behavior is like the 
"Zeno" or "watchdog" effect; there is no change at all from the initial state as iV — > oo. 

3.3.2 Storage with imperfect correction 

In any realistic situation we will have 5 > 0. For any finite value of 5, there is a 
non-vanishing probability of introducing a non-correctable error for each application 
of E+D, so it seems intuitively reasonable that there should be an optimum value of 
N. Also, for the case of storage, it makes sense to impose N5 < 1, or else the time 
taken for iV applications of E+D will exceed the time for which the qubit is stored. 
Practically, it is likely that N5 <C 1 would hold. 

The generalization of fllTD to iV equally spaced corrections is 

SNsc (n) = e ~ nN < TS (n e -(n-i) Kn ({i/N)-S)T _ ^ _ ^-^((i/at)-^ (24) 

For the simplest case of a decoherence rate always equal to K n , equating the derivative 
(with respect to N) to zero, keeping only the leading terms and rearranging to give 
the optimum N yields 

Note that this is consistent with our result (|22|), if we identify T opt with T/N. 

3.3.3 Transmission with imperfect correction 

Since the time for iV applications of E+D does not eat into T, but adds to it for 
transmission, there is not the absolute requirement that 5 < (1/N). However, for 
practical cases it is likely that N5 1, the same as for storage. The generalization of 
(|19|) to iV equally spaced corrections is 

t Nsc ( n ) = e ~ nN < TS (ne-^ KnT / N - (n - i)e~^ T / N ) N . (26) 
The optimum N is again given by (^), although there are differences at the next order. 
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3.3.4 iV-step optimization 



As expected, there is an optimum number N of corrections to apply when these pro- 
cedures themselves are imperfect. When NS <C 1, the optimum is the same for the 
storage and the transmission scenarios. 



It is interesting to substitute back the optimum iV of Q25D, to obtain the maxi- 
mum achievable transmission and storage success probabilities. At the first order of 
approximation they are equal and given by 



max(sNsc{n)) ~ max(tAr sc (n)) « exp —nn n T (2(n — l)5n n T) 



1/2" 

1 -riK n T(2(n-l)5K n T) l/2 . (27) 



Thus, within the simple framework used here, the minimum probability for a qubit 
state to incur an error in a total (storage or transmission) time T is approximately 
T(25(n-l)K n T) 1/2 . For cases of practical interest, where this probability is small 
(and so the success probabilities are close to unity), this is a good approximation. 



It is worth noting that, because of the square root in fl25|) , the optimum iV does 
not grow too quickly. For example, with phase noise at k — 10~ 5 (i.e. n = 3 and 
k 3 = 2k), and with A = 10 and a total time of T = 10 4 , the minimum qubit error 
probability (calculated using (p7|)) of 0.017 arises from applying N ~ 14 corrections. 
With isotropic noise at the same k (i.e. n = 5 and k 5 = 4k), and with A = 20 and an 
elapsed time of T = 10 3 , the minimum error probability of 0.016 is obtained with just 
two corrections. 



3.4 Errors and quantum jumps 

In all of this analysis we have been explicitly assuming that the influence of the envi- 
ronment produces errors of type (|j). It might be asked what the relation is between 
the general form of single-qubit errors given in (H) and the Lindblad master equation 
(§). At first glance they seem to have no resemblance to each other. The former is an 
abrupt, instantaneous change of state which occurs at random times; the latter is a 
continuous, deterministic equation for the density operator p. This is the more correct 
description of the system's evolution. In what circumstances can we approximate it by 

(§? 

If we represent the right-hand side of equation (|3|) by a superoperator C, then the 
master equation becomes 

|p = £p, (28) 
and given the density operator p(0) at some initial time we can formally solve for it a 
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time T later: 

p{T) = e CT p(0). (29) 
We can expand the right-hand side of this equation to get 



p(T) = e-^ T p(0)e^« T 

+ y [ T dt [e-^( r -*^ i e-^« fft p(0)e^ t Lte^^- t )" 
j Jo I 1 \ 

- k Jo Jo L 

ce^^te^^O^t^CT-t)] + . . . ? (3()) 



where 



H cS = H- l -J2L]L, (31) 

^ 3 



is a non-Hermitian "effective Hamiltonian" pl| , |22| , p3 



Already in this expansion we can see the relationship between the stochastic model 
of errors (|j) and the continuous master equation (0) Each term in ( |30D looks like a 
collection of instantaneous "jumps" (or errors) interrupting a continuous state vector 
evolution. 

However, it should be noted that this continuous evolution is not necessarily 
the desired evolution; the non-Hermitian component may produce unwanted effects, 
depending on the Lindblad operators. If we choose dephasing noise, so that L = 
i/k(T 2 , then the effective Hamiltonian is H c g = H — i(«/2)i, resulting merely in a 
renormalization of the state; in the case of spontaneous emission, however, we have 
L = and H c g = H — z(k/2)6" + <t_, which changes the relative weight of the |0) 

and |1) states. 

This sort of continuous error does not fit the error correction paradigm, and 
therefore cannot be completely corrected. However, all is not lost: if kT is very small, 
then it is possible to expand 

e - iA «* T « e" ii}T (l + KTO), (32) 

where O is a function of the commutator of H and L, and to first order in kT the error 
correcting algorithm will still work. 



Similarly, we see that the second and higher terms in ( pQ|) correspond to more 
than one "error" occurring during the time T; hence, error correction techniques for 
single errors will be ineffective for these terms. But again, for small kT, these terms 
will be of higher order, so correction is still beneficial. 
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The situation is somewhat more involved than this, however. For most models of 
a quantum computer the Hamiltonian will be time-varying. Thus, the time-evolution 
operators will not be the simple exponentials written in (j30|), but some more compli- 
cated operators. In the simple case where a gate is effected merely by "turning on" 
some Hamiltonian for a set period of time and then "turning it off" again, the time 
evolution operator for the operation of n gates would be 

f (T, 0) = e~ i6ntn ■ ■ ■ e-^e - ** 1 * 1 . (33) 
The expansion would then become 

p(t)= f eff (r,o)p(o)ft ff (r,o) 

+ E [ dt [ f es(T, t^tsit, 0)p(0)Tj fr (t, 0)L]f^(T, t) 

+ EE I V dt f dt ' [reff(r ) t)L i f eff (t,t , )4reff(t , ,o)p(o) 

; ^ JO JO 1 



xf} s (t', 0)Llf^(t, t f )L]ft s (T, t)]+-.-, (34) 

where the effective time-evolution operator T e s includes the non-unitary effects of the 
environment, just as in (|30|) . (Note that it is possible for the Lindblad operators to be 
different during the operation of each gate.) 

Let us now use this expansion to analyze the effectiveness of error correction in 
the presence of noise. For simplicity we will examine the three-qubit error correction 
scheme for dephasing noise. 

As we see from figure 1 the 3-qubit encoding scheme can be effected by a sequence 
of 5 Hamiltonians given in section |2.3j: 



f E (E, 0) = e -i H ^ e - iH ^ e - iH3ts e- iH2t2 e- iHltl , (35) 

while the decoding scheme is effected by applying the same gates in the opposite order: 

f D (D, 0) = e- i6lt5 e~ i62U e~ iA3t3 e- i&4t2 e- i66tl . (36) 

In this case, E = D = t\ + £2 + £3 + £4 + £5- The sum of these times is the time 
E + D = A defined above. The effects of the environment are summarized by three 
Lindblad operators of the form a/k<t 2 , one for each qubit. The effective Hamiltonians 
then become 

H cS = H-^l. (37) 

Assume that the qubits evolve for a time T between error correction steps, defined as 
the transmission scenario in Sect. |3.2.2| . The procedure is then as follows: the initial 
qubit is encoded into three qubits in a time E, evolves for a time T undisturbed, and 
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is then decoded in time D. Assume further that between the encoding and decoding 
stages, the Hamiltonian Hq = for the three qubits (interaction picture). The density 
operator is then 



p(T + A) 



e -3«(T+A)A (0) 

+ E / T IfDLjTE^TiVM dt 

7~i Jo J 



+J2 \f D f E {/\/2,t)L 3 f E {t^)m 



3 r A / 2 ~ 



xTi(t,0)LtTj(A/2,t)Tj]d* 



+ E / lT D (A/2,t)4f D (t,0)f EP (0) 



3 /"A/2 f 



xf E f j D {t,0)L]f j D {A/2,t)}dt 




(38) 



The third and fourth terms of this expansion represent the possibility of an error 
occurring during the encoding or decoding phase. Such errors cannot necessarily be 
corrected, and represent a loss additional to that from the higher order terms in the 
expansion. The first two terms represent the possibilities of no errors or a single 
correctable error occurring. 

The longer the time T between error correcting steps, the larger the higher-order 
terms, while the shorter the time T, the higher the proportion S = A/T spent encoding 
and decoding. Hence, there should be an optimal time T opt between error corrections 
as a function of k and A, which minimizes the total error rate, consistent with the 
result of (0). For dephasing or isotropic noise, (p2"D will hold exactly to lowest order. 

4 Numerical simulation 

Since the dimension of the Hilbert space of n qubits is 2 n , the density operator p in 
the n-qubit case can be represented by a complex Hermitian 2 n x 2 n matrix. This puts 
severe constraints on the memory of the computer used to simulate these systems. For 
n = 3 and n = 5, a direct numerical solution of the master equation (^j is feasible on 
a workstation. Simulating the master equation for larger values of n requires a much 
larger computer (more memory in particular), because of the exponential growth of 
the Hilbert space dimension. The difficulty of simulating such small n systems on a 
classical computer illustrates how much power would exist in a real quantum machine, 
where the computation would actually run in the Hilbert space. 
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Our n = 3 results (see below) were obtained by a straightforward integration 
of the density matrix, using a fifth-order Runge-Kutta algorithm for the numerical 
integration. For our n = 5 calculations, we have used an alternative approach, the 
quantum state diffusion method. This involves an average over sequential evolutions 
of a 2 n -dimensional quantum state, and therefore needs less computer memory. A 
workstation can probably handle up to about a dozen qubits if they are simulated 
this way. The drawback is that, because of the sequential runs required to construct 
good statistics, such a simulation may require a lot of CPU time. The n = 5 case 
can be handled by direct integration of the master equation and we have checked 
the accuracy of our quantum state calculations using direct integration of the master 
equation. Since the required computer memory grows like 2 2n for a direct solution of 
the master equation, but only like 2 n for a quantum trajectory simulation, the latter 
method can be used for values of n where the former would be impractical. 

Before we present our numerical results, we give a short description of quantum 
trajectories as numerical methods, with particular reference to quantum state diffusion. 



4.1 Quantum trajectory simulations 

The storage problem due to large density matrices can be overcome by unraveling the 
density operator evolution into quantum trajectories [19|, ^(], ^T], |22|, |23|, ^4|. Since 
quantum trajectories represent the system as a state vector rather than a density 
operator, they often have a numerical advantage over solving the master equation 
directly, even though one has to average over many quantum trajectories to recover 
the solution of the master equation. A single quantum trajectory can also give an 
excellent, albeit qualitative, picture of a single experimental run. 



We see from section |3J] that we can justify the use of the stochastic error models 
in sections 2 and 3, in spite of the continuous, deterministic character of the master 
equation itself. This type of treatment, in which the evolution of the density operator 
is written as a sum over many different stochastic evolutions of single wavefunctions, is 
called an unraveling of the master equation, and a single realization of these evolutions 
is a quantum trajectory. The unraveling of section 3.4 is often used in simulating 
quantum optical systems, where it is known as the "Quantum Jumps" or "Monte 
Carlo Wavefunction" approach |?T], [Z3] ]. 

The evolution of a single quantum jumps trajectory is given by the (Ito) stochastic 
differential equation 

m = _iH \^)dt - ~ £ (L)L 3 - (L)L^) \f)dt 
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+e(^^-iW^' (39) 

where the dNj are real stochastic differential variables which are except at cer- 
tain random times when they assume the value 1. These are independent, such that 
dNjdNk = 6jkdNj, and have a mean rate of jumps M(dNj) = (LjLj)^dt. Angular 
brackets denote the quantum expectation (G)^ = (ip\G\ip) of the operator G in the 
state The evolution between jumps is continuous and differentiable. The density 
operator is given by the mean over the projectors onto the quantum states of the en- 
semble. If the pure states of the ensemble satisfy the equation (|39|), then the density 
operator given by 

p(t) = M\m)(mi (40) 

satisfies the master equation (Q). From this it is clear that the expectation value of an 
operator O is given by 

Tr{6p} = M(ip\6\ip) . (41) 



Quantum jumps is a useful conceptual picture, but it is not the only unraveling 
of the master equation. It is convenient that we can use whatever unraveling we 
choose based solely on calculational convenience, as they are all equivalent to the 
master equation. Among the most important is the quantum state diffusion (QSD) 
equation of Gisin and Percival P0| . We have applied both jump and QSD equations 
to the problems considered in this paper. It turns out that to obtain good statistics, a 
significantly smaller number of trajectories need be summed when the QSD equation 
was used. We thus limit further discussion to the QSD equation, a nonlinear stochastic 
differential equation for a normalized state vector \ip): 

+ J2(L j -{L j ) i) ) IVR,-. (42) 

3 

The first sum in this equation represents the deterministic drift of the state vector 
due to the environment, and the second sum the random fluctuations. The d^j are 
independent complex differential Gaussian random variables satisfying the conditions 

Md£j = Md£id£j = , Md£*d£j = 5 i:j dt , (43) 

where M denotes the ensemble mean. A QSD trajectory is continuous, but not differ- 



entiable. If the pure states of the ensemble satisfy the QSD equation fl42|) , then the 
density operator given by (^0|) again satisfies the master equation (|3|). 

To simulate the QSD equation, we use a publicly available C+- 1- software library 
written by two of the authors [p5f| . The software uses object-oriented programming 
concepts to allow great flexibility in defining operators and states in Hilbert spaces 
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with arbitrary numbers of degree of freedom. As an illustration, we show how the list 
of Hamiltonians that defines the network effecting the encoding is implemented: 



const int n0fGates=10; 

Operator H[nOf Gates] = { A1+A2+A3, B023, C023, A4, B04, AO, 

B03+B14, B01+B24, A0+A4, B03 }; 

The operators in the list are implemented as follows (again we give as an example the 
operators A2 and B04 only): 

IdentityOperator id; 
SigmaX sx2(2) ; 
SigmaZ szO(O) ; 
SigmaZ sz2(2) ; 
SigmaZ sz4(4) ; 

Operator pr4 = . 5* (id+sz4) ; 
Operator prO = . 5* (id+szO) ; 

Operator A2 = (M_PI/2)*(sqrt(0.5)*(sx2-sz2) + id); 
Operator B04 = (M_PI) *pr0*pr4; 



4.2 Numerical Results 

The results obtained with the QSD method have been checked against those obtained 
by a direct integration of the master equation. The disagreement, of the order of 1-2% 
is purely statistical and is due to the finite number of trajectories used to build the 
average (around 200 trajectories for the three bit code and 200-400 for the five bit 
code). 

The simulations confirm the analytical results discussed in the previous section, 
both for the three-bit and the five-bit code. One measure of the efficiency of a quantum 
error correcting code is the mismatch between the decohered, corrected ensemble and 
the initial state, as defined in (0). This mismatch indicates how faithfully the initial 
state has been preserved in the face of noise. 

The mismatch m nec for a single qubit undergoing decoherence is defined by ( |i~4D 
for phase noise and ( |i~5l) for isotropic noise. Fig. |3] shows the isotropic noise mismatch 
of a single qubit. This is the benchmark to evaluate the efficiency of the five bit error 
correction code. 

A similar figure can be obtained by plotting the mismatch m ec of a qubit that has 
been encoded and later decoded. Instead of looking at the mismatch of a single qubit 
in contact with an isotropic noise reservoir for a time T, one encodes the qubit into 
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Figure 3: Mismatch between the initial and final states of a qubit in contact with an 
isotropic noise reservoir. The time scale starts at t = 20 to match the encoding and decoding 
times in comparing this figure with the numerical results for the case of a qubit with error 



correction. This follows the storage scenario of Sect 3.2.1 




Figure 4: Mismatch between the initial and final states of a qubit with a five bit error 
correction code. The time scale starts at t = 20 units to include the encoding and decoding 
times (A = 10 + 10 = 20). Each step (one or several gates) is effected in unit time, as 
indicated in Fig. |2[ 
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m aaaljtic m ec 




Figure 5: Difference between m ec and m ana i y tic, where m ec is the mismatch occuring in the 
numerical simulation and m ana i y ti C is the estimated mismatch from section 3. In both cases 
a single qubit is encoded into five and allowed to evolve in contact with a noisy environment, 
with noise occurring also during the encoding and the decoding phases. Time on the axis 
represents encoding and decoding time (i.e., A = 10 + 10 = 20 units of time) plus the time 
the encoded qubit is left interacting with the environment. 



five qubits (10 units of time), allows the five qubits to interact with the reservoir for 
T — 20 units of time, and finally decodes and corrects the qubit (10 units of time). The 
resulting reduced density operator is used to compute m ec via ([L3D . This mismatch, 
obtained by numerical simulation, is illustrated in Fig. f|. 

Alternatively, we can use the results of Sect. ^| and estimate the mismatch by 

m ana lytic = -(1 - s ac (5)), (44) 

where s sc (n) is defined in (|T7|). The agreement between this simple analytical model 
and the numerical simulations is very good, as illustrated in Fig. |5|. The maximum 
discrepancy is of the order of 10 percent for our range of parameters. 



Consistent with the analysis developed in Sect. |3.2.3| , one can identify in these 
numerical simulations a region of the k-T plane for which error correction is likely to 
help, despite noise occurring in the encoding and decoding stages. This can be clearly 
seen by looking (in analogy with Eq. |21]) at the positive values of log(m nec /m ec ). 

Another possible measure of the benefit of error correction is the difference m ec — 
m nec . Positive values of the difference indicate that error correction is worthwhile. 
However, this measure is of little use when both m ec and m nec go to zero, as their 
difference also vanishes. The log of the ratio does not have this drawback, and is 
therefore preferable as an indicator of where error correction is beneficial. The region 
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Figure 6: Surface and contour plots of log(m nec /m ec ), where m ec is generated numerically. 
The shaded area indicates values of (re, T) for which error correction is useful. This figure 
should be compared to Fig. 0. The thick line represents, for each given k, the optimum time 
between the E and D stages (cf. Sect. |3.2.3 , Eq. 22). 



of positive values of the log is represented by the shaded area in Fig. One notices, 
as expected, that for small enough k and large T, error correction is desirable (since 
\og(m nec /m ec ) > •<=>- m ec < m nec ). For comparison, Fig. ^ shows the same quantity 
where the analytical expression m ana i y u c of ( fHj) has been used instead of m ec . 



An exactly similar set of calculations can be done for three bit codes in the case 
of delocalizing noise, and similar behavior was observed. In both the three bit and five 
bit cases, there is a section of the k-T plane where error correction remained beneficial 
even in the presence of noise during encoding and decoding; and for low values of 
the environmental interaction strength k, there was an optimal time between error 
correction steps. This is consistent with the result obtained by Chuang and Yamamoto 

m 



5 Conclusions 



From both the analytical arguments and the numerical simulations, we see that error 
correction can prove worthwhile even in the presence of noise during encoding and 
decoding. For a given strength of the environmental coupling, there is an optimal rate 
at which error correction should be performed, and for a given time of storage there is 
an optimal number of error correction steps. 
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Him nec lm analytic) 

b g( m ase /m analytic ) . 




Figure 7: Same as Fig. |, but here we plot log(m nec /m ana i ytic ) instead of log(m nec /m ec ). 
The discrepancy between the numerical calculations and the analytical arguments is illus- 
trated in Fig. ||. 

The simulations presented in this paper are among the first to treat both the 
execution of gates and the influence of the environment realistically, in the sense that 
the operation of gates takes a finite amount of time during which noise continues to act 
on the system. Moreover, models of the noise were used which correspond to common 
environmental effects in atomic and optical physics. 

While the theory of error correction has moved rapidly, it is unlikely that circuits 
involving many qubits will be experimentally realized soon. Systems of a few qubits 
thus remain of great interest. Three-bit and five-bit error correction are among the 
first circuits that might be experimentally implemented, and hence our results should 
be of relevance to near-future experiments in this field. 

We also have seen that quantum trajectories provide a practical technique for 
simulating systems with multiple qubits. This may prove particularly useful in treating 
systems with many qubits, where solving the full master equation is impractical due 
to the large size of the Hilbert space. 

These simulations could be improved and extended in many ways. The Hamil- 
tonians used to represent the gates were chosen for convenience rather than reflecting 
any particular physical system. It would be useful to get closer to the actual physics 
of proposed quantum computers, such as the linear ion trap of Cirac and Zoller []26| . 
In the same way, the coupling to the environment might differ for different gates. It 
would be straightforward to include these effects. 
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There are other interesting problems in quantum computation which might be 
studied by techniques like those of this paper. The recently proposed fault-tolerant 
error correction schemes are far more complicated than the ones treated in this paper. 
They would be beyond the reach of direct numerical simulation of the master equa- 
tion with present computers, but may well prove amenable to a quantum trajectory 
approach. 

Quantum computation still faces many hurdles before becoming reality. But it is 
far too early to say that the ingenuity of those working in the field is not sufficient to 
overcome them. 
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